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ABSTRACT 

We have measured the three time delays for the quadruple gravitational lens RX J1131 — 1231 using 
two seasons of monitoring data. The short delays between cusp images are A-B = 11.98j^[ p days 
and A-C = 9.61 1 \ ®l days. The long A-D delay for the counter image is not as precisely determined 
because of the season gaps, but the data suggest a delay of — 87±8 days. The short delays are difficult 
to explain using standard isothermal halo models of the lensing potential, which instead prefer A-B 
and A-C delays of ~ 1 day for reasonable values of the Hubble constant. Matching the cusp delays is 
possible by adding a significant (~ 5 x 10 10 M Q ) amount of matter nearly coincident (~ 0'.'05 South- 
East) with the A image. Adding such a satellite also helps improve the quasar and lens astrometry 
of the model, reduces the velocity dispersion of the main lens and shifts it closer to the Fundamental 
Plane. This is suggestive of a satellite galaxy to the primary lens, but its expected luminosity and 
proximity to both image A and the system's bright Einstein ring make visual identification impossible 
even with the existing HST data. We also find evidence for significant structure along the line of sight 
toward the lens. Archival Chandra observations show two nearby regions of extended X-ray emission, 
each with bolometric X-ray luminosities of 2-3 x 10 ergs s _1 . The brighter region is located A6 rj 
153" from the lens and centered on a z — 0.1 foreground cD galaxy, and the fainter and presumably 
more distant region is 4-5 times closer (in angular separation) to the lens and likely corresponds to 
the weaker of two galaxy red sequences (which includes the lens galaxy) previously detected at optical 
wavelengths. 

Subject headings: gravitational lensing: individual (RX J1131— 1231) 



1. INTRODUCTION 

Gravitationally lensed quasars are useful tools for 
studying the matter distribution of intermediate red- 
shift (0.1 < z < 1.0) galaxies. Robust estimates of 
the total (luminous plus dark) lensing mass usually fol- 
lows from just measuring the image positions. How- 
ever, quantifying the presence of halo substructure, con- 
straining the mass profile's radial distribution, or pars- 
ing the relative contributions from stars and dark mat- 
ter requires additional information. For example, the 
flux ratios of lensed images are sensitive to the pres- 
ence of halo substructure (e.g., iMetcalf fc Zhaol 120021) 
and the incidence of so-called "anomalous" ratios at ra- 
dio wavelengths is consistent with CDM predictions that 
a few percent of the halo mass should be left in satellites 
(Dalal ~ Kochanekl 12002ft . At optical wavelengths, in- 
terpreting such anomalies in light of stellar microlensing 

1 Based on observations obtained with the Small and Moderate 
Aperture Research Telescope System (SMARTS) 1.3m operated 
by the SMARTS Consortium, the NASA/ESA Hubble Space 
Telescope, and the Spitzer Space Telescope. HST observations 
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is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with HST program 9744. The Spitzer 
telescope is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under a contract with NASA. Support 
for this work was provided by NASA through an award issued 
by JPL/Caltech. These observations are associated with Spitzer 
program 20451. 

2 Department of Astronomy, Ohio State University, 
Columbus, OH 43210; nmorgan@astronomy.ohio-state.edu, 
ckochanek@astronomy.ohio-state.edu, xinyu@astronomy.ohio- 
state.edu 

3 Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138; falco@cfa.harvard.edu 



implies that < 25% of the surface mass density at the 
image locations must be in a smooth (non-stellar) f orm 
<)Schechter fc Wambsganall2002l iKeeton et al.ll2006|h 

Another key observation is time delays. The time 
delays between lensed images measure a combination 
of the Hubble constant and the surface mass den- 
sity (n) = (£) /S c of the lens at the radius of the 
lensed images, w ith H cx (1 — (k)) to lowest order 
ijKoch anck 2002). For th e six simpl e time delay lenses 
that we have at p resent, [K ochanek fl2002fl noted that 
four (PG1115+080 . ISchechter et alill 9971 iBarkanalll 997 



llmpev et al.1 



1998; SBS1 520+5 30, 



Burud et al.fp 002b; 



B1600+434, IBurud et all 1200(1 iKoopmans ; et al J I200C : 
and HE2149-2745. IBurud et aLll2002al) yield H„ predic- 
tions consistent with the HST Key P roject result of 72 
km s _1 Mpc -1 ijFreedman et al.l|2001[) only if the lensing 
galaxies have constant M/L profiles (that is, declining 
rotation curves at the system's Einst ein radius). The re- 
maining two s ystems (HE1 104-1805 lOfek fc Maozl2003l 
HE0435-1223. lKochanek et al.ll2u0q) have measured de- 
lays consistent with the Key Project only with mass pro- 
files that drop off slower than isothermal (that is, rising 
rotation curves at the locations of the lensed images). In 
the case of HE0435-1223, the best-fit model yields an av- 
erage mass convergence at the system's Einstein radius 
that is 20% larger ((k) ~ 0.6 ± 0.05) than expected for 
an isothermal halo. Taken at face value, the range in (k) 
estimates imply a heterogeneous mix of mass profiles for 
lensing galaxies, with the likely explanation that we are 
seeing the effects of the lensing environment. Lens galax- 
ies located at the center of their respective groups will 
naturally have an extended dark matter halo (yielding 
flat or rising rotation curves), while lens galaxies located 
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in the field or group periphery will have mass distribu- 
tions dominated by their stellar components (falling ro- 
tation curves). Under this interpretation, it becomes im- 
portant to not only expand the sample of systems with 
measured time-delays in order to sample as wide a range 
of environments as possible, but also to characterize the 
lensing environment for systems with measured delays as 
thoroughly as possible. 

In this paper, we present our time-delay measurements 
and explore the lensing environment for the quadruple 
lens RX J1131— 1231 (hereafter J1131). Our original mo- 
tivation was to use the measured delays and assumed H D 
to estimate the lens surface mass density at the J 1131 
image positions. Coupled with the robust total mass es- 
timate obtained from the image separations and the mea- 
sured luminosity distribution of the lens obtained from 
flST/NICMOS data, this would allow us to estimate the 
radial mass profile of the lens and separate the luminous 
and dark contributions to th e total lensing mass (a s done 
recently for HE 0435-1223; [Kochanek et all l2006). Un- 
fortunately, the measurement of unexpected delays for 
this system makes such an approach unfeasible. Instead, 
we simply focus on constructing a plausible macromodcl 
that can reproduce the delays while assuming an isother- 
mal halo for the main lens. 

This system was originally identified as a ROSAT X- 
ray source and su bsequently reveale d to be a quadruply- 
imaged quasar bv tSluse et a il pOO^ . The lensing geom- 
etry consists of a background z s = 0.658 quasar lensed 
by a foreground zi — 0.295 elliptical galaxy, and there is 
a prominent E instein ring evid ent even in ground-based 
optical data. iSluse et alJ l)2003|) also detected variability 
on the order of 0.3 magnitudes from the total integrated 
system flux over an eight month baseline. Overall, the 
system's demonstrated variability, wide image separation 
(A to D separation of 3'.'2), prominent lens galaxy and 
quad plus ring morphology made J1131 an attractive tar- 
get for long-term monitoring and followup studies. In |3 
and <J3] we describe our monitoring program and time- 
delay measurements and comment on the microlensing 
variability observed in the system over the course of the 
campaign. High-resolution HST/ ACS and NICMOS im- 
ages of the lens and immediate environment are described 
in 21 an d evidence for at least two clusters along the 
line of sight are presented using far IR (Q and X-ray 
(S0 observations obtained with the Spitzer and Chandra 
Space Telescopes. We then explore several iterations of 
parametric and non-parametric lens models and discuss 
their implications in In summary, we find it diffi- 
cult to simultaneously explain the system geometry and 
measured time delays using a single profile for the main 
lens galaxy, but find it possible if a significant amount 
of substructure is present near one of the cusp images. 
Finally, we summarize our results and describe avenues 
for future work on this system in SJH1 

2. LENS MONITORING 

The monitoring data for J1131 were obtained as part 
of the ongoing program described by iKochanek et alJ 
(2006); details of the reduction pipeline and lens fitting 
procedures can also be found in that reference. Briefly, 
we monitor approximately 25 lensed quasars at cadences 
of 1-2 times per week, with typical observations of three 
2-5 minute exposures per epoch in either Sloan r or John- 



son R. Each lens is assigned up to 5 nearby reference 
stars which serve as local flux calibrators and as shape 
templates for the analytic PSF model. Point sources 
are modeled as a combination of three elliptical Gaus- 
sians and extended objects (lens galaxies) are modeled 
using Gaussian approximations to de Vaucouleurs pro- 
files (which allows for rapid analytic convolutions with 
the seeing disc). The relative astrometry for each system 
is fixed using observations from the CASTLES 4 HSTpra- 
gram, and relative photometry for the system images are 
obtained by minimizing x 2 residuals between the data 
and analytic model. 

The monitoring data for J1131 were obtained using the 
queue-scheduled SMARTS/CTIO 1.3m telescope in the 
Southern Hemisphere and the MDM 2.4m, WIYN 3.5m, 
and APO 3.5m telescopes in the Northern Hemisphere. 
The bulk (94%) of data were taken with the SMARTS 
1.3 m using the dual-be am optical/IR ANDICAM cam- 
era (pePov et al.ll2003j) and covering the 2003-2004 and 
2004-2005 seasons for a total of 101 epochs. These obser- 
vations consisted of three 5 minute i?-band exposures and 
six 2.5 minute J-band exposures at each epoch. The IR 
data are considerably less sensitive than the optical data 
and we focus only on the i?-band light curves here. Ad- 
ditional observations include three Sloan r-band epochs 
consisting of 1.5-2 minute exposures using the SPIcam 
CCD on the APO 3.5m, two i?-band epochs of 2-3 minute 
exposures using the WTTM camera on the WIYN 3.5m, 
and two i?-band epochs of 3 minute exposures using the 
8K Mosaic CCD on the MDM 2.4m. The APO, WIYN 
and MDM light curves were calibrated by interpolating 
the quasar and reference star data points onto the re- 
spective SMARTS curves, yielding star offsets of —0.01, 
0.00 and —0.02 mags and quasar offsets of 0.08, —0.02 
and 0.04 mags, respectively. We only kept data with see- 
ing better than l'/7 FWHM, or about half the maximum 
image separation. 

FigureQ]shows the four quasar light curves and Tabled 
gives the quasar magnitudes and reference star differen- 
tial variability at each epoch. Over the two seasons, each 
image varied by as much as 0.7 mag. Significant corre- 
lated variability on timescales of weeks is also apparent 
for the three brighter images. 

3. TIME DELAYS MEASUREMENTS 

The light curves are modeled as a sum of two Legendre 
polynomial series. The first series, of order N src , models 
the intrinsic source variability of the quasar and is the 
same for all four quasars images up to a time delay shift 
for each pair. The second series, of order N^, is differ- 
ent for each image and models the different macrolens- 
ing magnifications (0 th order) and the uncorrelated mi- 
crolensing variability for each quasar image (1 st order 
and above). Exact prescriptions for the fitting forms can 
be found in IKochanek et aTl l)2006tL The merit function 
for each light curve is minimized by differentiating with 
respect to the polynomial coefficients and solving the re- 
sulting linear set of equations, which yields estimates for 
the time delays between each image, a lightcurve for the 
intrinsic source variability, and the relative microlensing 
lightcurves for three of the four images. 

4 Due to the career moves of the investigators, CASTLES is no 
longer an acronym but rather just the name of the survey. 
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Fig. 1. — Light curves for images A-D. The B and D curves are offset for clarity. The best fit models are for N STC = 30, Nu = 4 for each 
separate season for images A, B and C, and for -/V src = 40, Np = 3 for the joint seasons for image D. Differential photometry for 5 nearby 
reference stars is shown in the bottom panel. 



The A, B and C images show relative delays of about 
1-2 weeks, much shorter than the length of each observ- 
ing season. Since all the signal when solving for the ABC 
delays will come from same-season data, we fit the two 
season's ABC light curves simultaneously (that is, we 
solve for a single set of time delays constrained by both 
seasons) but model each season's lightcurve separately. 
(We discuss modeling the D curve below). The F-test is 
used to decide when further increases in N src and iV M do 
not yield improvements to the fits. We find that most 
of the source-curve improvements occur up to iV src = 10, 
beyond which the higher source-curve complexity does 
not significantly improve \ 2 ■ Significant improvement in 
the microlensing model occurs from — (macromodel 
magnifications only and no microlensing) to = 4, 
arguing for microlensing variations that are somewhat 
more complex than just a simple linear trend (which 
would be described by = 1). We adopt A src = 30, 



= 4 as our standard model for the A, B and C images, 
where the larger source curve order should lead to conser- 
vative uncertainties in the delay estimates. This model 
gives A-B, A-C and B-C delays of 11.98±i;p 9.61+^ 
and — 2.20l}'g4 days (ler), respectively. As a check on 
the results, we also fitted each season's light curve as 
a stand-alone dataset using a separate code that uses 
brute-force (Powell) minimization, and we find delays 
consistent with the quoted error bars. Figure [2] shows 
the ABC light curves after shifting by the relative delays 
and subtracting the microlensing and macromodel mag- 
nifications, leaving just the source variability observed in 
each image. 

Fitting for image D's light curve is less precise. In 
addition to the lower signal-to-n oise, D is expecte d to lag 
the other images by ~ 100 days (Slu se et al.l2003() . which 
would shift about half of D's data points into the season 
gaps of the other images. This leaves less overlapping 
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Fig. 2. — Overlapping curves for images A (open squares), B 
(triangles) and C (filled squares) after shifting by the best-fit time 
delay and removing the macromodel magnification and microlens- 
ing variability from images B and C. The much noisier D curve is 
not shown. Errorbars are suppressed for clarity. 



data to measure a delay, but we can improve the available 
signal by modeling the two seasons as one continuous 
light curve. Fixing the relative ABC delays at the values 
determined above, we use a N SIC = 40, = 3 model to 
obtain an A-D delay of — 87±8 days (1<t). The reported 
best-fit is driven by aligning the broad valley in image 
D at MJD w 3125 days with the similar valley in the 
ABC curves at MJD « 3040, with most of the season 
two drop in D's brightness shifted into the ABC season 
gap. If this feature is not correctly matched, then the 
absence of other distinguishing features in the D curve 
would mean that the A-D delay is unconstrained with 
the extant data. If the feature is correctly matched, then 
the D delay is measured to an accuracy of about 10%. 

Figure shows the residual microlensing light curves 
after subtracting the appropriately time-delay shifted 
source variability (arbitrarily set to be image A's light 
curve) from each quasar image. Both B and C show 
a drop in brightness relative to A toward the end of the 
first season, and image B shows a rapid rise in brightness 
at the end of the second season. Image D does not show 
significant variability apart from a gradual (Am ~ 0.05 
mag) decline during both seasons, although the lower 
signal to noise complicates the analysis. 

It is likely that image A is responsible for the first 
season drop seen in the B and C curves since we only 
measure the differential effects. This would imply a net 
increase in A's brightness by ^0.2 mag due to microlens- 
ing toward the end of the first season. The differences 
between the B and C microlensing light curves are evi- 
dence for a similar amount of microlensing in one or both 
of these images as well. For example, B's upswing at the 
end of season two is not seen in image C. Also, image 
B spends most of season two ~0.2 mag fainter than its 
season one average, while image C maintains more or less 
the same average flux from one season to the next. 

One succinct way to quantify the effe cts of microlensing 
is to look at the lensing cusp relation ijMao fc Schneider! 



1998), which states that the sum of the signed image 
magnifications is zero for the triplet of images formed 
when a source is near a cusp caustic (images A, B 
and C for J1131). The relation can be phrased in 
terms of the observed image fluxes (Fa, Fb, Fc) by di- 
viding through by the unknown source flux, yielding 



\F A + F B + F C \/{\F A \ + \F B \ + \F C \). While 



the cusp relation for a source asymptotically close to a 



cusp has i? c 



0, in practice i? C usp as defined will 



increase by several tenths as the source is placed far- 
ther from the cusp, not to mention the a dditional per- 
turbing effects of micro- and millilensing l|Keeton et alJ 
2003, hereafter KGP). For a system configuration such as 
J 1131, one can place a strong upper limit on the baseline 
-Rcusp value (before micro- or millilensing perturbations, 
that is, assuming a perfectly smooth potential model) 
of <0.1 (KGP). This value is significa ntly lower than the 
R cusp w 0.35±0.03 measured from the lSluse et all (2003) 
V-band discovery data, and lead KGP to classify J1131 
as a "cusp-anomaly" system. 

Anomalous i? CU sp values indicate the presence of sub- 
structure in the lens galaxy, but it can not discern 
whether the anomaly is due to microlensing by stars or 
due to millilensing by substructure (either dark matter 
clumps or luminous satellites). However, only microlens- 
ing can change i? C usp on observable timescales, so track- 
ing changes in i? C usp tracks the microlensing contribu- 
tions. The bottom panel of Figure shows i? C usp com- 
puted using the time-delay corrected light curves (shift- 
ing to the image A time frame). i? C usp is ~0.23 during 
most of the first season but begins to drop, becoming 
less "anomalous" with respect to the KGP estimate to- 
ward the season gap. At the same time, B and C drop 
in brightness and the i? CU sp value remains at ^0.14 for 
the bulk of the second season. The res olved Chandra X- 
ray observations of J1131 reported bv iBlackburne et all 
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Fig. 3. — (Top): The microlensing light curves for 
RX J1131 — 1231. Shown are the light curves for images A-D after 
subtracting the time-delay shifted A model curve as a model for 
the source variability. (Bottom): R C usp for images A, B and C 
formed using the best-fit curves in Figure after shifting by the 
respective time delays. 
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Fig. 4. — PyDrizzled HST/ ACS /-band image of RX J1131— 1231 and the surrounding field. Prominent nearby galaxies are labeled. 
Concentric circles mark the 10" and 20" radii centered on the lens galaxy. The inset shows the inner 10" region at the same scale but at 
10 X higher contrast. 



( 2005) have R cusp = 0.75 ± 0.05 on MJD = 3108 due to 
the significantly dimmed flux of component A (see their 
Figure la). Taken together, the SMARTS and Chandra 
data suggest that the saddlepoint A image was in a de- 
magnified microlensing state for much of season one and 
varied closer to the unperturbed macromodel magnifica- 
tion heading into season two. 

The HST observations of J1131 described in the next 
section yield respective R cusp values in the V, I and 
H filters of 0.26 ± 0.06, 0.17 ± 0.05 and 0.08 ± 0.01. 
Assuming the baseline i? C us P value of < 0.1 (KGP), 
the trend with wavelength agrees with the notion that 
microlensing preferentially affects smaller emission 
regions of the quasar accretion disk (under the assump- 
tion that longer wavelengths correspond to lower disk 
temperatures, which then correspond to larger radii 
from the central black hole). However, it is important 
to remember that the cusp anomaly may be due to 
the presence of substructure on scales comparable to 
the image separations. For example, models with a 
satellite galaxy near one of the cusp images (explored 
in H7.3(l predict an i? C us P value of about 0.5 even before 
any microlensing effects. In such a case, the offsets 
from -Rcusp J; 0.1 cannot be interpreted as strictly a 
sign of microlensing, but the variations in R cusp seen 
in time and across wavelength still point to abundant 
microlensing for J 1131. 



4. HST OBSERVATIONS 

High-resolution images of J1131 w ere obtained with 
the ACS/Wide Field Camera fWFC: iFoTTTTaTI ITflfl^ 
and NICMOS detectors on board the HST as part of 
the CASTLES imaging program (PID 9744). Both of 
th ese public data sets ha ve also been recently analyzed 
bv lClaeskens et all l|2006fl , and we make comparisons be- 
tween the separate reductions below. The optical ACS 
images were taken with the F814W (/-band) and F555W 
(V-band) filters on 22 and 24 June 2004, respectively. 
Exposures were 33 minutes through each filter using a 
5-point dither pattern. The data were reduced using the 
IRAF/CALACS package as part of the "on-the-fly" re- 
processing at the time of download. Subsequent cosmic- 
ray rejection, geometric correction, and image combina- 
tions were performed using the standard PyRAF pro- 
grams available for ACS data reduction. The infrared 
NICMOS images were taken with the NIC2 camera and 
F160W (iJ-band) filter on 17 November 2003, and con- 
sisted of two sets of four-point dither MULTIACCUM 
exposures for a total on-source integration of about 89 
minutes. The NICMOS data were reduced using the NI- 
CRED software described in lLehar et al.1 l)2000jf . 

Figure 0] shows a 1' field of view surrounding J 1131 
from the stacked /-band data. The environment within 
10" of the lens is devoid of galaxies, and the nearest 
galaxies to the lens (G3, G4 and G7) are 15-20" distant. 
The Figure 0] inset shows the inner 10" region centered 
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Fig. 5. — Combined ffi?T images of RX J1131 — 1231. (a): ACS/V-band image, (b): Same as (a), but after subtracting the best-fit point 
source, lens galaxy and ring model, (c): NICMOS/-ff-band image, (d): Same as (c), but after subtracting the best-fit point source, lens 
galaxy and ring model. Contrast in both residual panels stretches from — 7% to +7% of the respective peak lens galaxy flux. 



on Jl 131 at 10 x higher contrast to emphasize the lack 
of nearby structure. The nearest galaxies with compara- 
ble magnitudes to the lens are Gl and G5, roughly 30" 
distant toward the North-North- West. 

In Figures and [SJ: we show the V- and iZ-band 
closeup images of J 1131. A complex pattern of lensed 
knots from the quasar host galaxy is present, presum- 
ably because the relatively low redshift host galaxy is 
a star-forming spiral galaxy. There is also a faint fifth 
ob ject, component E in Figure El (denoted component X 
in Iblaeskens eFal][2006), roughly 0'.'5 North of the lens 
galaxy. The object shows a clear Airy ring in the iZ-band 
data, so it is unresolved. 

To model the light distribution, we first tried fitting 



the 7Z-band data using a seven component model: five 
point-sources for images A-E, a de Vaucouleurs profile 
for the lens galaxy, and a lensed de Vaucouleurs profile 
for the quasar host galaxy to model the continuous por- 
tion of the Einstein ring emission. For the fitting, we 
used the IMFIT pa ckage written by B. McLeod (see, e.g., 
iLehar et alJ 120 00). The lensing potential was modeled 
as a singular isothermal sphere plus external shear us- 
ing just the iZ-band image positions as constraints for 
the purpose of mapping the host galaxy onto the im- 
age plane. We also tried exponential and Gaussian pro- 
files for the quasar host, but the de Vaucouleurs profile 
yielded the lowest residuals. Figure EJi shows the resid- 
uals after fitting a circularly symmetric de Vaucouleurs 
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profile for the lens galaxy and reveals the lens to be elon- 
gated roughly along the A-D line. Allowing an elliptical 
de Vaucouleurs profile for the lens does well at modeling 
the galaxy ellipticity (Figure^) but still leaves extended 
residuals near the center of the lens with a peak flux of 
about 12% of the unsubtracted galaxy peak flux. To 
model this extra emission, we added an exponential pro- 
file centered on the de Vaucouleurs model. This addition 
did a better job at accounting for the bulk of the ex- 
tra emission, although there are systematic residuals in 
the core of the lens galaxy (~10% of the unsubtracted 
galaxy flux; Figure EJ) that are suggestive of a slightly 
misaligned center for the de Vaucouleurs and exponential 
profiles. 

The inability to mo del the main le n s with a single 
profile was noted by IClacskcn s et alJ l)2006[) as well. 
The two-component galaxy model naturally suggests a 
disk-f bulge morphology, except that we find the de Vau- 
couleurs model to be a better fit to the elongated struc- 
ture than to the bulge. For an SO galaxy, one would 
expect the opposite, namely an exponential disk and a 
de Vaucouleurs bulge. We tried to coax the models to 
assign the exponential profile to the disk by using judi- 
cious choices for the initial conditions, but the models 
consistently preferred a de Vaucouleurs disk and expo- 
nential bulge. The de Vaucouleurs component of the 
galaxy model is best characterized by an effective radius 
R e = 1"58±0'.'29, axis ratio of 0.51±0.02, and major axis 
orientation of —62° ± 2°. The exponential component of 
the galaxy has a best-fit scale length of 0'.'20 ± 0.01, axis 
ratio of 0.79 ± 0.04, and major axis orientation 33° ± 
4°. 

For the V- and /-band data, we fixed the astrometric 
and structural components of the system at the values 
from the ii-band eight component model and solved for 
the relative photometry. Table lists the results for all 
three filters. The relative positions for the A, B and 
C images agree withi n the quoted errorbars with the 
iClaeskens et al.l l|2006) positions, but we do find a small 
offset of 0'.'009 (3cr) for the D image and a much larger 
offset of 0'.'029 (15a) for the main galaxy position. The 
lens offset between the two reductions is about half an 
ACS pixel and may originate with the different light pro- 
files used when modeling the main lens galaxy (de Vau- 
couleu rs + Exponential for our fit, single Sersic profile 
for the IClaeskens et alJl2006l fit) . 

The colors of the system components are plotted in 
Figure El Components A and B are three-quarters of a 
magnitude redder in I — H compared to C and D, but 
all have nearly identical V — I colors. Component E is 
significantly redder by about 1.5 mag in V — I than the 
other four quasar images. The reddening vector assum- 
ing a Ry = 3.1 Galactic extinction law is given by the 
arrow in Figure EI and shows that E is unlikely a red- 
dened copy of a fifth quasar image. Its I — H color is too 
red by several tenths of a magnitude for a late-K/early- 
M dwarf as well. Its color is consi s tent w ith the primary 
lens galaxy, and Clae skens et alJ 1)200 6) suggested that 
E may be an unresolved satellite galaxy. 

The main lens galaxy has total (exponential + de Vau- 
couleurs) colors in V — I and I — H of 1.84 ±0.11 and 
1.83 ± 0.20, respectively. These are bluer by ~0.2 mag- 
nitudes than expected for az/ = 2.5 burst galaxy model 
observed atzw 0.3 (see Figure EJ, but the difference is 
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Fig. 6. — Colors for RX J1131-1231 system components A- 
E and lens galaxy components Gd (de Vaucouleurs component) 
and Ge (exponential component). Color tracks for the dwarf main 
sequence (heavy line) and a zj = 2.5 burst galaxy model (thin 
line) are also drawn. The arrow gives the reddening magnitude 
and direction assuming a Ry = 3.1 G alactic extinction law and 
wavelength-dependent parameters from Cardel li et alJ J198SI1 . 



not unreasonable compared to the spread in galaxy col- 
ors observed in oth er lensed systems at similar redshifts 
l)Rusin et al.ll20 03) . As for the exponential contribution 
to the galaxy model, its magnitude is fainter than the de 
Vaucouleurs component by 2.1 mags, 1.5 mags and 1.2 
mags in the V, I and H filters, respectively. Its colors are 
redder than the de Vaucouleurs component, as expected 
for a bulge-like component to the lens galaxy. However, 
as already noted, the morphology is not as expected for 
a true bulge and the spatial offset between the "bulge" 
and "disk" components evident in Figure [S^ is puzzling. 

5. SPITZER OBSERVATIONS 

We also observed J1131 with the Infrared Array Cam- 
era (IRAC) on board the Spitzer Space Telescope on 11 
June 2005. Observations consisted of a 36-point dither 
pattern for each of the four IRAC bandpasses (3.6, 4.5, 
5.8, and 8 fim), covering an area of roughly 7f 5 x 7.'5 
centered on the lens. The data received from the Spitzer 
Science Center had passed through the standard IRAC 
pipeline (consisting of flat-field corrections, dark subtrac- 
tion, and linearity and flux calibrations) and the resulting 
Basic Calibrated Data (BCD) files served as our input for 
MOPEX drizzling. The BCD images were combined us- 
ing a final scale of 0'.'863 pixel -1 , and Figure [7\ shows the 
[3.6] image of J 1131 and surrounding field. 

Even after drizzling onto the finer pixel grid, the J1131 
A-D separation spanned only four pixels. The lack of 
spatial resolution complicates any analysis of the lens 
components and we instead focused on the IR colors of 
field galaxies. After identifying exten ded sources in the 
field using the SExtractor package ijBertin &: Arnoutsl 
1996L ver. 2.2.2), we searched for red sequences (e.g., 
Glad ders fc Yeel I2000D in the IR color-magnitude dia- 
grams indicative of overdense regions of elliptical and 
SO galaxies (groups or clusters). The resulting CMDs 
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were suggestive of a low-redshift (z < 0.3) over density of 
galaxies with [3.6]-[4.5] colors of ~ 0.1, but the concentra- 
tion was not defined enough to estimate a redshift. Since 
the Sp itzer observations were obtained, I Williams et alJ 
(2005) have indeed identified two red sequences for the 
Jl 131 field using deep VRI observations obtained with 
the KPNO Mayall 4m and CTIO Blanco 4m telescopes 
(see their Figure 8). The dominant red sequence has a 
centroid (denoted by the eastern bullseye in Figure 0) 
located 2.'1 East of' J1131 G and «45" to the South- 
South- West from a bright elliptical galaxy identified in 
the mid- infrared images (GO in Figure 01. Galaxy GO 
is suggestive of a cD g alaxy associated with the primary 
IWilliams etafl l)2005|) red sequence, implying significant 
structure along the line of sight toward Jl 131 (addi- 
tional evidence supporting this claim is presented using 
the archival X-ray observations describe d below). The 
second , weaker red sequence detected bv IWilliams et al.l 
( 2005) has a flux-weighted centroid 25" North- West of 
J 1131 G (western bullseye in Figure [7J and includes the 
lens galaxy itself. 

6. CHANDRA OBSERVATIONS 

Given the evidence for structure along the line of sight, 
we re-analyz ed the archival 10 ks C handra observation 
of J1131 (see iBlackbur neet alJ l2005) taken on 12 April 
2004 using the ACIS-S3 detector, with the goal of iden- 
tifying X-ray emissio n from halos associated with the 
IWilliams et all pjbH) red sequences. The smoothed X- 
ray emission (after excising the X-ray flux associated 
with the lensed quasar) is shown as contours on top of 
the IRAC [3.6] image in Figure [7| We identified two ex- 
tended sources of X-rays in the field aside from the lens 
flux analyzed by iBlackburne et all l)2005l) . The first is 
the bright emission centered 153" North-East of J1131 
(at ll h 32 m l?4, -12° 31' 6" J2000; centroid accuracy 
of approximately 2"). The X-ray position is coincident 
with the large elliptical galaxy (GO) identified in the mid- 
infrared images and within an a rcminute of the red se - 
quence flux centroid identified bv IWilliams et all |2005). 
This galaxy (and two others within three arcminutes) 
has a redshift measured from the Las Campanas Redshift 
Survey of z = 0.10 and th e region was identifi ed as an 
optically-selected group bv iTucker et al.1 l)200flK Thus, 
the red sequence, the cD galaxy and the extended X-ray 
emission confirm the presence of a foreground cluster at 
z = 0.1. 

We extracted the A CIS spectrum o f the extended emis- 
sion and used XSPEC <|A rnaudl!996|) to fit it with a ther- 
mal plasma model modified by Galactic absorption and 
obtained a good fit to the data with a x 2 = 22 for 41 de- 
grees of freedom. We estimated an unabsorbed 0.4-8 keV 
flux of (7.5 ± 0.8) x 10~ 13 ergs cm" 2 s" 1 , a bolometric 
X-ray luminosity of (3.1 ± 0.5) x 10 43 ergs s _1 assuming 
the cluster is at z — 0.1, and a temperature of 1.13±0.03 
keV. The temperature and luminosity are roughly con- 
sistent with typical L-T relations fe.g.. | W u et al]|1999f 
[ Helsdon fc Ponmanll200rl IXue fc WulEo Ort IRosati et"aTl 
2002) for groups and clusters, providing further confir- 
mation of the cluster redshift. 

The second extended source of X-ray emission is South- 
West of the lens, but it is fainter and the analysis is 
complicated by other point sources in the field. To pro- 
ceed, we did a general search for faint extended emis- 




FlG. 7. — Smoothed X-ray contours plotted on top of the IRAC 
Channel 1 (3.6 /ira) image. X-ray contours are logarithmic ex- 
cept for the first two contour levels, which mark the 2cr and 3<r 
confidence levels. The eastern (western) bullseye marks the raw 
ce ntroid of the primary (secondary) galaxy red sequence detected 
bv lWilliams et alJ 120051) . 

sion f ollowing the procedures used in iDai & Kochanekl 
()2005|) . We used the CI AO tool wavdetect to identify 
point sources in the ACIS-S3 field and then replaced 
them with nearby background regions to create a soft 
band (0.5-2 keV) image. We smoothed this image with a 
Gaussian of width a = 25" to find a source 33" Southwest 
of the lens (ll h 31 m 50?l, -12° 32' 23", J2000). Because 
of the low (~ 3a) significance of the detection and the ef- 
fects of masking the overlapping emission from the lens, 
the position uncertainty is comparable to the smoothing 
scale. Since this flux overlaps the lens and is close to the 
centroid of the weaker red galaxy sequence detected by 
Williams et al. 2005 (western bullseye in Figure [TJ), we 
tentatively identify the emission as from the lens group 
at redshift of z; = 0.295. The measured count rate in the 
0.5-2 keV band is 0.004 cts s" 1 within a circle of radius 
25", corresponding to an unabsorbed flux of 1.7 x 10~ 14 
ergs cm~ 2 s _1 . If we use a /3 model to extrapolate the X- 
ray flux to a larger aperture and assume a temperature of 
T ~ 1.5 keV, then we estimate that the bolometric lumi- 
nosity of this second cluster is 2 x 10 43 ergs s -1 . Thus, 
the halo associated with the lens group is comparable 
in bolometric luminosity to the halo associated with the 
low redshift group, but would require 10 x the integration 
time to characterize at the same signal-to-noise. 

7. LENS MODELS AND INTERPRETATIONS 

In this section, we describe our mass modeling ef- 
forts aimed at reproducing the HST geometry and the 
SMARTS time delays. We discuss three model varia- 
tions. First, we look at both isothermal and variable- 
slope power-law halo models for the main lens galaxy 
( H7.1|l . but in general find that such models cannot ex- 
plain the long (~ 10 days) A-B and A-C time delays. 
Next, we try a non-parametric mass model of the system 
( H7.2[l . In general, such models add no information be- 
yond that available in parametric models, but the asym- 
metric mass reconstruction suggests a satellite galaxy or 
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other substructure near image A may be required to ex- 
plain the anomalous delays. Finally, we return to para- 
metric models fi J7.3[) consisting of the primary isother- 
mal halo for the main lens and a secondary halo for the 
presumed substructure near image A. Such models can 
plausibly reproduce both the long A delays and the image 
astrometry at the expense of consuming the remaining 
degrees of freedom. 

7.1. Basic Isothermal Models 

We first tried to model J1131 using the standard sin- 
gular isothermal ellip soid embedded in an external shear 
l|Kassiola fc Kovnerl Il993t lKorman n"eraT] Il99l . For 
constraints, we used the HST/ ACS -ff-band positions 
(Table |5J with error bars of 3 mas for the quasar po- 
sitions and the three SMARTS time delays. We looked 
at three model sequences: model SISx uses a spherical 
halo centered on the lens galaxy plus an external shear; 
model SIEx allows for halo ellipticity plus shear; and 
model SIEx-l- allows the halo center to move with respect 
to the galaxy center but constrained with errorbars of 5 
mas. The astrometry and time delays give 13 constraints, 
so the models have 6, 4 and 4 degrees of freedom, respec- 
tively. We ignored the flux ratios. All models were min- 
i mized in t he im age-plane using the gravlens software 
of lKeetonl (j2001jl and adopting a (ft m ,n A ) = (0.3,0.7), 
H Q = 70 km s -1 Mpc" 1 cosmology. The model results 
are summarized in Table [3] 

None of the models gives a statistically acceptable fit 
to the observations; the best-fit SIEx+ model has a fi- 
nal x 2 /dof of 74.1/4. The simpler SISx and SIEx models 
cannot reproduce the image astrometry, with rms values 
between the observed and predicted quasar locations of 
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Fig. 8. — Schematic of RXJ1131— 1231 environment. Labels X- 
ray 1 and 2 denote the large and small X-ray emissions detected in 
the Chandra observations, and Red Seq. 1 and 2 denot e the two red 
sequence centroids identified by Williams ct al. (2005); the dashed 
lines marks the la spread in the external shear direction for the 
SIEx+ model; the pentagons mark the locations of neighboring 
galaxies identified in the HST/ ACS data. The diameter of the X- 
ray 2 circle is the X-ray smoothing scale of 25". Symbol sizes for 
the galaxies are proportional to the square root of their ACS I- 
band fluxes. The lens galaxy /-band flux (not depicted by symbol 
size) is nearly twice that of Gl. 



33 and 13 mas, respectively (final x 2 /dof values of 538/6 
and 138/4, with more than half of the \ 2 budget com- 
ing from the quasar positions). The SIEx+ model fits 
the quasar positions essentially perfectly (image rms of 2 
mas) at the expense of shifting the halo center by 12 mas 
(4<r) from the measured lumino us galaxy positio n . The 
situation is worse if we use the iClaeskens et all <|2006|) 
image and galaxy positions, which yield x 2 /dof values of 
534/6, 158/4, and 130/4 for the SISx, SIEx and SIEx+ 
models, respectively. A substanti al external shear is also 
required for both the Table l2*land lClaeskens et ahl 1)200611 
positions - greater than 10% for all models - and points 
at —83° East of North for the SIEx+ case (where the 
convention is for the shear direction to point towards or 
away from the perturber.) There is no obvious perturber 
along the shear direction (see Figure |HJ) ; a galaxy with 
roughly the same velocity dispersion as the lens would 
produce a 10% shear at 9", but the nearest galaxies of 
significance (G3 and G4) are 20" away and at the wrong 
position angle. Galaxy G8 lies at the correct position 
angle, but is much too distant to produce a 10% shear 
given its brightnesses and that of the main lens galaxy 
(see Figure 0}. The shear does not point toward either 
of the X-ray halos detected in the Chandra data, or to- 
ward either of the re d sequence centroids identified by 
Willi ams et alJ l)2005|) . This likely indicates that either 
there is no single external perturber dominating the shear 
or that the external shear is attempting to compensate 
for an inadequacy in the main lens model. 

The more troublesome problem is that all models fail to 
reproduce the measured time delays. The SIEx+ profile 
predicts A-B and A-C delays of 0.98 and 1.23 days (as- 
suming H a = 70 km s -1 Mpc -1 ) compared to the mea- 
sured values of roughly 12 and 10 days. Formally these 
are 9a and 5er discrepancies. Note that the predicted B- 
C delay of 0.25 days for the SIEx+ model agrees within 
2<T with the measured value of — 2.20t 1 'g4 days, so the 
problem is chiefly with the A image. The situation can- 
not be corrected by simply rescaling H a , since matching 
the A-B and A-C delays would require an unrealistic 
Hubble constant of w 10 km s -1 Mpc -1 . 

The long A delay is difficult to explain. There are many 
examples of model perturbations leading to systematic 
uncer tainties in time delay predictions (see iSchechteil 
120051 for a review), but these typically alter the de- 
lay predictions at the 10% level and not the factor of 
10 n eeded here. For exam ple, the mass-sheet degener- 
acy (jOorenstein et alJll988f ) rescales the image delays by 
At oc (1 — (k)) by superimposing a mass-sheet of surface 
density (n). One might expect a convergence of (k) ~ 
0.1-0.3 from matter associated with the nearby X-ray 
halos (e.g., as with syste ms Q0957+561. iBarkana et all 
Il999t RXJ 0911+0554, iBurud et alJ Il998|) . but "this 
would uniformly shorten all delays by 10%-30%, and not 
preferentially lengthen the A delay. To a similar degree, 
small variations in the radial exponent of the halo density 
profile ( ~ 10%) alter dela y predictions by a comparable 
fraction <|Kochanekl 120021. again of insufficient magni- 
tude to explain the A delay. An example of an extreme 
change in the density slope is to use a de Vaucouleurs 
constant mass-to-light ratio model with an effective ra- 
dius matching the value measured from the NICMOS 
data. Such a model can match the astrometry as well as 
the SIEx+ model provided the de Vaucouleurs center is 
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Fig. 9. — PixeLens ensemble-average mass map of 
RX J1131— 1231, as constrained by the system astrometry and 
ABC time delays. The mass contours are logarithmically spaced at 
factors of 2.5 in surface density, beginning with the outer contour 
at k ps 0.16. 



allowed to float, but the predicted cusp delays are still 
~ 1 day. 

Abandoning the assumption of isothermality, we also 
explored a range of elliptical power-law profiles with 
three-dimensional density profiles falling as shallow as 
1/r to as steep as 1/r 4 . We found that highly ellipti- 
cal profiles steeper than 1/r 3 could produce long cusp 
delays provided the halo was orientated along the A-G 
line. Because of the high ellipticity, such models concen- 
trated more mass near A than near B and C, increasing 
the Shapiro contribution to the A delay. However, the 
required ellipticity to produce an A-B delay of even half 
the observed value was flatter than e = 0.9. Such a mass 
distribution would indicate a prominent disk galaxy, and 
this is clearly inconsistent with the morphology of the 
lens observed in the NICMOS data. 

7.2. Non-Parametric Models 

Lacking a successful macromodel, we turned to non- 
parametric models us ing the PixeLens software of 
iSaha fc Williams! l)2004() . As noted above, such models 
do not add any more information than present in para- 
metric models since they are grossly underconstrained, 
but they may indicate where the parametric models are 
qualitatively wrong. PixeLens works by reconstructing a 
pixelated mass map of the lens using the image positions 
and time delays as constraints. For J 1131, we recon- 
structed both symmetric and asymmetric mass models 
over a 6" x 6" grid centered on the lens with a "pixrad" 
of 11, yielding around 530 separate mass elements. Since 
PixeLens assumes negligible errors for all input values, 
we tried models both with and without the more uncer- 
tain D delay but found little qualitative change in the 
final mass map. 

The reconstruction for the asymmetric mass distribu- 
tion is shown in FigureEl (For the symmetric reconstruc- 
tion, the eastern half of the structure is mirror-imaged 
on the western half, producing an unrealistic four-leaf 



clover mass distribution). There is a clear protrusion 
surrounding image A. PixeLens's solution for increasing 
the A delay is to have significant variation in the angu- 
lar distribution of matter along the three cusp images 
by increasing the surface mass density around image A 
(ft ~ 0.9) as compared to the surface mass densities at 
images B and C (n ~ 0.1-0.2). Qualitatively, this is sim- 
ilar to our result found for a highly elliptical disk-like 
halo overlapping the A image. However, the asymmetric 
reconstruction stresses that what is chiefly needed is a 
higher surface mass density at A compared to the other 
two cusp images, suggesting that a satellite galaxy or 
other undetected substructure at this location may be 
responsible for the anomalous delay. 

7.3. Substructure Near Image A 

There is no visual sign of a satellite galaxy near im- 
age A in the HST images or the corresponding subtrac- 
tion residuals. A red satellite galaxy would be most ev- 
ident in the F160W image, but the significant Einstein 
ring emission and PSF residuals complicates any visual 
search near the lensed images. To explore if a perturber 
could reproduce the cusp delays, we added an SIS mass 
model over a grid of positions surrounding image A and 
looked at how each model adjusted for the J1131 time 
delays. The main lens galaxy was again modeled using 
the SIEx+ profile, and we fitted both the astrometric and 
delay constraints as before. We found that a perturber 
with an Einstein radius of ~ 0'.'2 and placed within ~ O'.'l 
South-East of the A image could lengthen the A-B and 
A-C delays by several days, but such models consistently 
formed two additional images of the background quasar. 
The perturber's proximity to image A means that its crit- 
ical curve always maps into a source-plane caustic near 
the source position, creating a six-image region inside 
the primary lens astroid caustic. While one of these ex- 
tra images is always highly demagnified and located at 
the perturber position, the other typically forms at twice 
the perturber-A separation with a magnification ~ 10% 
that of image A. Such an image would be easily identified 
in the HST data. 

We can avoid the formation of extra images by adding 
a core radius s to the perturber's profile, such that the 
projected density is given by 



(e.g., IKassiola fc Kovnerl Il993l IKormann et~aH I1994J) . 

The exact criteria needed to prevent extra images depend 
on the perturber's core size and mass scale. Consider a 
model where the main lens (modeled as an SIE) and per- 
turber (modeled as a softened isothermal sphere, or IS) 
are colinear along the x-axis with the main lens at the 
origin and the perturber at a distance x = x' . Including 
an external shear term (with shear angle B 1 measured 
from the x-axis), the potential 4>(x) is 

2 

4>{x) = fam(x) + 4>is(x - x') + ^- cos(20 7 ). (2) 

The appearance of extra images first becomes possible 
(that is, a critical line and corresponding six-image caus- 
tic are first formed) when a solution to the lens equation 
exists at x = x'. This translates into a condition on the 
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Fig. 10. — Optimal perturber positions for the best-fit sub- 
critical (ft = 1.4s, ft = 1.6s, b = 1.8s), critical (ft = 2.0s) and super- 
critical (ft = 2.1s, ft = 2.2s) models. The three contours outline the 
Ax 2 = 2.6, 4.3, 6.2 (lcr, 2<r, and 3cr for two parameters) contour 
parameters for the critical case. The location of image A is marked. 
For scale, the dotted box shows the size of one ACS pixel centered 
on image A. 



potential of 



dx 21 



or, after substitution, 



2s 



1. 



(3) 



(4) 



where we have neglected a shear term that is always much 
smaller than unity for reasonable models. No extra im- 
ages are possible when b' /2s < 1 and multiple images 
are possible (provided the source is inside the six-image 
caustic) otherwise. 

We looked at a range of perturber models: two super- 
critical (b — 2.2s and b = 2.1s), a critical (b = 2.0s) and 
three sub-critical (b — 1.8s, b — 1.6s and b — 1.4s). Dur- 
ing the minimization, a large penalty term was added to 
the merit function whenever the source position crossed 
into the six-region caustic, insuring that only four-image 
models were accepted. We also constrained the orienta- 
tion of the main lens halo using the measured NICMOS 
position and associated error, which helped stabilize the 
ellipticity of the primary lens. This approach gave a final 
tally of 12 parameters (two for the source position, seven 
for the main lens, and three for the perturber) and 14 
constraints for two degrees of freedom. Table summa- 
rizes the result for the critical perturber. 

Figure I1UI shows the optimal perturber positions for 
the model sequence and the confidence contours for the 
critical (b — 2.0s) case. For all models, the perturber po- 
sition was tightly constrained to a region between 0703 
to O'.'ll South-East of image A. In Figure ITTI we show 
the resulting y^/dof for each model as the perturber's 
core radius is increased from s = 0" (the singular case) to 
s = 0'.'3 (about 1.3 kpc at the lens redshift). Note that for 
the two super-critical models, \ 2 rapidly increases once a 
threshold core radius is crossed, with roughly 90% of the 



increase coming from the quasar positions. The abrupt 
increase is due to the growing size of the six-image caus- 
tic as the perturber's core radius and potential depth are 
increased, which displaces the source position from its 
optimal location inside the six-image region. This means 
a super-critical perturber is highly unlikely since it re- 
quires a fine-tuned source position in order to match the 
quasar astrometry while remaining outside the six-image 
region. The critical and sub-critical models avoid the 
formation of the six-image caustic and consequentially 
show no abrupt increase in \ 2 - 

As seen from Figure^] the critical model is marginally 
preferred over the sub-critical cases. There are two no- 
table improvements for the optimal critical model com- 
pared to the single halo model considered in H7.ll First, 
both the quasars and galaxy astrometry are now fitted 
essentially perfectly: x 2 /dof — 0.3/2 compared to 13/4 
for the SIEx+ case, with virtually all (99%) of the x 2 
budget now coming from the time delays. Second, the 
cusp delays are lengthened considerably: A-B, A-C and 
A-D delays of 5.90, 7.88 and —80 days, respectively. The 
remaining model parameters are quite reasonable as well, 
with a nearly round main lens galaxy (ellipticity of 0.14) 
and small external shear (4%). Thus, it is possible to 
improve both the overall astrometry and the time delay 
predictions by adding a perturber halo nearly coincident 
with image A, provided the halo is soft enough to avoid 
the formation of extra images. 

The A-B delay is still a factor of two smaller than mea- 
sured (a 4<7 discrepancy). Since the Shapiro contribution 
to the time delays depends on differences in effective po- 
tential between images, one simple way to fine-tune the 
delays is to vary the radial slope of the perturber. With 
a variable slope a, the projected surface density n(r) can 
be written as 



1 



U2-a 



2 ( s 2 +r 2)l-a/2- 



(5) 



where again s is the core radius and the normalization 




perturber core radius s (") 

Fig. 11. — Goodness of fit for two super-critical (ft = 2.2s and 
ft = 2.1s), critical (ft = 2.0s) and sub-critical (6 = 1.8s) perturbcrs 
near image A. Minimums for the ft = 1.6s and ft = 1.4s cases are 
also noted by the two tick marks. 
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FlG. 12. — Potential strength (ft') and halo slope (a) for a non- 
isothermal perturber. The dotted line marks the potential strength 
for the just-critical perturber, using each model's core radius and 
best-fit halo slope. 



is chosen such that b is the system's Einstein radius in 
the limit of an isothermal profile (a — > 1) and vanishing 
core radius. For small core radii, the effective potential is 
approximately given by (p(r) ss b' 2 ~ a r a /a 2 . Generating 
longer A-B and A-C delays requires a larger difference in 
the A-B and A-C effective potentials, which is possible 
for steeper than isothermal (a < 1) mass profiles. 

To test if a non-isothermal slope could more closely 
match the delays, we allowed the potential depth b 1 and 
halo slope a to vary independently while again looking 
at a range of core radii s. This modification consumes 
the remaining two degrees of freedom and we expect a 
perfect fit for a plausible model. Indeed, we found that a 
halo slope with a = 0.7-0.8 could reproduce both the as- 
trometry and measured time delays essentially perfectly 
although with a clear degeneracy between the size of the 
core radius and the halo slope. For the s = Of! 2 model 
(summarized in TableEJ, the predicted delays were 12.0, 
13.8 and -85 days for A-B, A-C and A-D, respectively, 
with an overall x 2 /dof of 4.7/0. It is interesting that 
we still find x 2 > even for zero degrees of freedom. 
The only x 2 contribution for all s > 0'.'15 models is from 
the A-C delay, although it is more instructive to look at 
the discrepancy using the B-C difference. The predicted 
B-C delays (between +1.0 and +1.8 days for the range 
of core radii considered) is in the opposite sense of the 
measured value (—2.20^64 days). The delay ordering 
amon g images is a model-ind ependent feature of a lens 
(e.g., iSaha fc Williams! 120081) . so if the sense of the de- 
lay is measured incorrectly then even a model with zero 
degrees of freedom will be unable to reproduce the in- 
correct delay sense. All models considered in this paper 
predict B trailing C (B-C > 0), implying that images 
A and B are the "merging" image pair normally seen in 
inclined-quad configurations and that the overall image 
ordering is CBAD. Therefore, if the measured C delay is 
off by 2-3 a such that the measured B-C delay becomes 
positive, then it might explain the non-zero x 2 - 



Figure lT2l shows the degeneracy between the perturber 
slope a and core radii s, as all models in the figure yield 
X 2 /dof values between 3/0 and 5/0. We also plot the 
best-fit potential depth b' for the perturber and the con- 
dition for the just-critical case (dotted line), which for 
the power-law profile is 



2\s 



2-a 



1. 



(6) 



The total x 2 /dof gradually drops from 4.7/0 at s = 0'/15 
to 3.0/0 at s = 0'.'45 with again the only significant con- 
tribution coming from the C delay. Over this range, b' 
if fairly constant, ranging from 0'.'25 to 0'.'28, and the 
halo slope increases from a = 0.7 to 0.8. Note that 
for s > 0'.'15, the perturber no longer sits on the just- 
critical line as seen for the optimal (b — 2s) isothermal 
model. Thus a steeper-than-isothermal profile not only 
reproduces the long A delay but also naturally avoids the 
formation of extra images. 

7.4. Mass of the Perturber 

The cylindrical mass for the surface density given in 
Eq.|3]increases as M cy \(r) oc r a , so formally the total per- 
turber mass in the adopted model is unbounded. Adopt- 
ing a cut radius comparable to the A-C or A-B image 
separations (wl'/2, or 5.3 kpc at the lens redshift) gives 
an enclosed perturber mass of 4.7 x 10 10 M©. As a com- 
parison, the mass enclosed by the primary lens galaxy out 
to the same cylindrical radius (centered on its respective 
profile) is 4.4 x 10 11 M©, so the perturber has < 10% of 
the mass of the main lens at comparable radii. This is 
fairly significant and one might expect some visual evi- 
dence if it is indeed a satellite galaxy. We can estimate a 
plausible luminosity to the perturber by scaling its veloc- 
ity dispersion with that of the main lens and using the 
iFaber fe Jacksor] |l976) relationship. For an enclosed 
mass that scales with velocity dispersion as M (R) oc o 2 , 
the factor of 10 difference in mass translates into a factor 
of 3.2 in velocity dispersion, or a factor of 100 difference 
in luminosity assuming L « a 4 . Such a galaxy would 
be extremely difficult to see given the ring emission and 
light from image A, even after image subtraction. For ex- 
ample, scaling from the observed count-rate of the main 
lens measured in the NICMOS image, the predicted peak 
count rate from the perturber is only ~0.008 cnts s _1 . 
The rms count rate within a 0l!A radius around image 
A after subtracting the quasar images and ring emission 
(Figure |3f) is 0.058 cnts s _1 , which is larger than the 
peak signal expected from the satellite galaxy by a fac- 
tor of seven. Thus even if the perturber were a bona fide 
satellite galaxy, it is doubtful we would detect it in the 
existing HST images. 

7.5. Placement on the Fundamental Plane 

One check on the models is to compare the pre- 
dicted velocity dispersion of the main lens galaxy with 
the value expected from the Fundame ntal Plane (FP; 
iDiorgovski fc Davisll987llDressleill987l) . Almost all lens 
galaxies lie on the present-day FP provided one allows 
for passive luminosity evolution with redshift ijTreu et alJ 
120011 I2002t IRusin et al.l 12003^1 . When studying grav- 
itational lenses, it is convenient to phrase this evolu- 
tion in terms of the galaxy's B-band mass-to-light ratio 
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Fig. 13. — M/Lb evolution of the J1131 lens galaxy (arrows 
offset in redshift for clarity) from the present-day FP, computed 
using combinations of the total de Vaucouleurs plus Exponential 
galaxy magnitude, just the de Vaucouleurs component magnitude, 
and the two velocity dispersion estimates described in the text. 
T he other 18 p o ints a re evolutions obtained for lens galaxies listed 
in IRusin et al.l 120031) that have spectroscopi cally measur e d lens 
redshifts. The solid lines denote the best-fit Rusi n et a T] J2003) 
slope of Alog(Af/L) = —0.54 ± 0.04 using 28 gravitational lens 
galaxies. 



(M /L) ft measured i nside the system's Einstein radius 
fsee lTreu et al.ll2001l for a thorough discussion). Param- 
eterizing to first-order in redshift as Alog(Af/L)s = 
az, studies of grav itational len s gala xies find values 
between a » .5 ijRusin et alJ l2003|) and a « 0.7 
llTreu et alJ 120051). Simila r results a re found for cluster 
ll van Dokkum et al.lll998|) and field ijvan Dokkum et ahl 
1200 1|) galaxies out to z ~ 0.8 that are not associated with 
lensing. 

Using a, Zf = 2.5 starburst SED to calculate the neces- 
sary fc-corrections and evolution effects, the restframe B- 
band surface brightness for J1131 G is 20.62 ± 0.05 mag 
arcsec -2 using the total (de Vaucouleurs plus exponen- 
tial) galaxy magnitude and 20.83 ± 0.03 mag arcsec -2 
using just the de Vaucouleurs galaxy magnitude. The 
system's velocity dispersion is either 353 ± 2 km s -1 or 
337 ± 2 km s- 1 for the SIEx+ or SIEx+ plus a = 0.72 
perturber model. Both values correspond to extremely 
bright galaxies, 6L* and 5£» respectively for c* = 225 
km a -1 . Such galaxies are rare but not unheard of. The 
IBernardi et"aiT<|2003aj) SDSS sample of ~9000 field ellip- 
ticals contains galaxies with velocity dispersions as large 
as 400 km s" 1 . Although a a — 350 km s -1 galaxy is 200 
times less common than a a « 170 km s -1 galaxy found 
at the peak of the elliptical velocity dispersion function 
ijSheth et al.ll2003|K the cross-section for lensing scales 
as (7 so massive galaxies are preferentially favored. Fi- 
nally, the galaxy intermediate-axis effective radius mea- 
sured from the NICMOS data gives a physical size of 
R e = 5.0 ± 1.2 kpc. 

The evolution rates required to place J1131 G onto 
the present-day FP are plotted in Figure Four es- 
timates are computed, using the possible combinations 
of the two surface brightness estimates and the two ve- 



locity dispersion estimates. Also plotted are rates com- 
puted for 18 other gravitational le nses with m e asure d 
lens and source redshifts listed in Rusin et~aT] (J2003). 
Contrary to the IRusin et alJ l)2003|) sample of lens galax- 
ies, all J1131 G evolution rates are inconsistent with 
passive luminosity evolution, with only the halo plus 
perturber model and total (de Vaucouleurs plus expo- 
nential) galaxy magnitude rate consistent with no evo- 
lution at all at the la level (a = 0.02 ± 0.05). An- 
other way to view this is that the galaxy's inferred ve- 
locity dispersion is simply too large for its measured 
luminosity. Adopting M* = -19.79 + 51og(fe) ± 0.04 
llMadgwick et alJ l2002). 7 = 4.0 ± 0.25 l|Bernardi et alJ 
2003bfando7 = 225±22.5 km s" 1 , the Faber- Jackson re- 
lationship yields a = 274±28 km s -1 using the combined 
de Vaucouleurs plus exponential magnitude. This esti- 
mate is considerable lower than the value of a — 337 ± 2 
km s -1 obtained from the halo plus perturber model, 
suggesting that there may be an additional source of 
mass contributing convergence inside the Einstein ring 
besides the main lensing galaxy. Part of the solution 
likely lies in the mass convergence associated with the 
overlapping X-ray halos detected in the Chandra data 
and the unknown object component E, if it is truly a 
satellite galaxy to the primary lens. Both effects would 
lower the inferred value of the main lens velocity disper- 
sion by accounting for some of the inferred mass inside 
the system's Einstein ring. 

8. SUMMARY AND CONCLUSIONS 

We have measured the time delays for the quadruple 
lens RX J1131 — 1231 based on two seasons of monitor- 
ing data from the SMARTS/CTIO 1.3m telescope. The 
short delays for the cusp images are A-B = 
and A-C = 9.61^^57 and the long delay for the counter 
image is A-D = -87 ± 8. The A-B and A-C delays 
are an order of magnitude larger than expected using a 
standard elliptical model of the lens galaxy. One plau- 
sible explanation arises from a satellite galaxy or other 
dark substructure about 0'.'05 from the central cusp im- 
age, although some fine-tuning is required to match the 
observed properties of the system. In particular, for 
the power-law mass models explored here, the perturber 
must have a core radius of > 0'.'2 (> 1 kpc) to avoid the 
formation of extra images of the background quasar, a 
steeper than isothermal mass profile (a ~ 0.7; p oc r~ a ) 1 
and a significant total mass within the A-C and A-B im- 
age separations of ~ 4.7 x 10 10 M Q (about 10% of the 
main lens galaxy for comparable radii). Although mas- 
sive, such a satellite galaxy would have a luminosity only 
~1% of the main lens, making it undetectable in the ex- 
tant data given the perturber's proximity to the bright 
quasar A image and the prominent Einstein ring emis- 
sion. 

Ancillary evidence for such a perturber comes from the 
improved astrometry of t he quasar positions and primary 
lens galaxy. As noted bv lClaeskens et ail l)2006|) . using a 
standard isothermal halo for the main lens can reproduce 
the relative quasar positions but not th e position of the 
lensin g galaxy. This was addressed by iClaeskens et alJ 
(2006) by adding an m = 4 octupole term (character- 
istic of boxy/disky isodensity contours) to the lensing 
potential, although they admit that the physical nature 
of the term is difficult to interpret. Comparable devia- 



14 



Morgan et al. 



tions from pure elliptical models have been ruled out for 
four other l ens ga laxies using the shapes of Einstein rings 
l|Yoo et al.ll2005|) . Here, we have shown that including a 
perturber near image A allows for an essentially perfect 
fit to the quasars and galaxy astrometry while simulta- 
neously explaining the anomalous cusp delays. 

We have als o re-analyzed the C handr a X-ray data first 
presented by iBlackburne et al.1 l)2005j) and found evi- 
dence for two X-ray halos along the line of sight toward 
J1131. The brighter halo is coincident with a foreground 
[z = 0.1) cD galaxy 2.' 5 North-East of the lens and cor- 
responds to the stronger of two galaxy red seq uences 
detected in the optical by iWilliams et alJ 1)20(151) . The 
fainter halo is much closer to the lens and presumably 
includes J1131 itself, given the halo's proximity to the 
lens and proximity to t he weaker of the two galaxy red- 
sequences identified by IWilliams et al.1 l)2005|h 

Overall, the data suggest a very complicated lensing 
environment for J 1131. The lens galaxy is likely part of 
az« 0.3 group or cluster of galaxies with sufficient mass 
to support a ~ 1.5 keV X-ray halo; there is a foreground 
cluster of galaxies along the line of sight as evidenced by 
a second X-ray halo of similar temperature and centered 
on a prominent z = 0.1 cD galaxy; a significant mass 
perturbation near the A image is required to explain the 
anomalous time delays for the cusp images; modeling 
the surface brightness profile of the main lens requires 
at least two distinct profile shapes, suggestive of a more 
complicated primary lens than a simple de Vaucouleurs- 
shaped elliptical; and finally, there is a faint, point-like 
object 0'.'5 North of J1131 G with colors similar to those 
of the main lens galaxy but otherwise of an as-of-yet un- 
determined nature. These complications make it difficult 
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at present to use the measured time delays to either es- 
timate the Hubble constant or to constrain the overall 
mass profile of the main lens. The monitoring program 
has revealed extensive (~ 0.3-0.4 mag) microlensing vari- 
ability in the optical on time scales of months, so future 
microlensing studies of the system are promising. 

Progress with future mass modeling will at least re- 
quire measuring the central velocity dispersion of J1131 
G to break the mass-sheet degeneracy between the main 
lens galaxy and X-ray halos identified in the Chandra 
data. Such information will likely help to reconcile the 
lack of apparent M/L evolution for the main lens as well, 
since it will provide an estimate of the halo convergence 
inside the Einstein ring that is currently assigned to the 
lens galaxy. Deep imaging and spectral observations of 
component E are obviously required as well, but may 
prove challenging given its overall faintness, extremely 
red color, and proximity to the main lens. 
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TABLE 1 Lightcurves for RXJ 1131-1231 



HJD x 2 /dof Comp. A 



Comp. B 



Comp. C 



Comp. D 



ref. stars 



Telescope 



2988.811 
2995.804 
3002.819 
3013.819 
3021.790 
3028.765 
3035.815 
3039.758 
3042.783 
3046.734 
3050.769 
3052.918 
3055.713 
3058.707 
3063.712 
3066.739 
3067.760 
3070.686 
3074.674 
3079.825 
3081.715 
3084.852 
3088.616 
3092.574 
3103.735 
3104.662 
3107.612 
3111.596 
3114.620 
3129.597 
3132.609 
3133.718 
3136.473 
3139.503 



SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 

APO 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
EIGHTK 
SMARTS 

WTTM 
SMARTS 
SMARTS 
SMARTS 

WTTM 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
SMARTS 
EIGHTK 
SMARTS 
SMARTS 



1.29 
5.62 
10.58 
2.33 
6.08 
4.21 
3.93 
2.89 
0.76 
1.41 
7.66 
1.50 
2.08 
4.69 
4.85 
1.14 
3.05 
0.99 
1.25 
17.12 
3.84 
13.07 
3.85 
3.52 
2.04 
6.17 
4.61 
2.37 
3.75 
1.85 
0.77 
5.61 
1.36 
2.64 



2.379 ± 0.009 

2.465 ± 0.010 
2.383 ± 0.007 
2.398 ± 0.008 
2.495 ± 0.009 

2.460 ± 0.007 
2.454 ± 0.007 

2.461 ± 0.008 

2.479 ±0.012 
2.440 ± 0.009 
2.443 ± 0.007 
2.505 ±0.006 
2.491 ± 0.009 
2.457 ± 0.007 

2.466 ± 0.008 
2.518 ±0.011 

2.480 ± 0.008 
2.460 ±0.011 
2.471 ±0.010 
2.431 ± 0.005 

2.467 ± 0.007 
2.294 ± 0.004 
2.437 ± 0.007 
2.405 ± 0.007 
2.393 ±0.009 
2.333 ± 0.004 
2.378 ± 0.007 
2.383 ±0.008 
2.375 ± 0.007 
2.300 ± 0.009 
2.315 ±0.011 
2.329 ±0.011 
2.321 ±0.010 
2.266 ±0.007 



2.233 ± 0.009 
2.220 ± 0.009 

2.321 ± 0.007 
2.388 ± 0.008 
2.342 ± 0.009 
2.369 ± 0.007 
2.385 ± 0.007 
2.404 ± 0.008 
2.334 ±0.011 
2.361 ± 0.009 
2.361 ± 0.006 
2.365 ± 0.006 
2.328 ± 0.008 
2.324 ± 0.007 
2.377 ± 0.007 

2.322 ±0.010 
2.350 ± 0.008 
2.337 ±0.010 
2.311 ±0.010 
2.237 ± 0.005 
2.263 ± 0.007 
2.158 ±0.004 
2.254 ± 0.007 
2.275 ± 0.007 
2.222 ± 0.008 
2.168 ±0.004 
2.241 ± 0.007 
2.244 ± 0.008 
2.211 ±0.007 
2.160 ±0.008 
2.100 ±0.011 
2.103 ±0.011 
2.086 ± 0.009 
2.099 ± 0.007 



3.192 ±0.014 
3.088 ±0.013 
3.147 ±0.009 

3.270 ±0.013 
3.144 ±0.013 

3.271 ±0.010 
3.284 ±0.010 
3.232 ±0.012 
3.264 ±0.021 
3.249 ±0.015 
3.216 ±0.009 
3.341 ± 0.008 
3.213 ±0.012 
3.271 ±0.010 
3.271 ±0.011 
3.243 ±0.015 
3.256 ±0.012 
3.232 ±0.018 
3.241 ±0.017 
3.225 ±0.007 
3.230 ±0.010 
3.168 ±0.005 
3.224 ±0.011 
3.192 ±0.010 
3.173 ±0.014 
3.177 ±0.005 
3.142 ±0.010 
3.171 ±0.012 
3.157 ±0.010 
3.117 ±0.014 
3.105 ±0.017 
3.065 ±0.015 
3.074 ± 0.014 
3.060 ±0.010 



4.560 ±0.037 
4.380 ± 0.027 

4.495 ± 0.021 
4.595 ± 0.034 
4.499 ± 0.028 

4.496 ± 0.021 
4.512 ±0.021 
4.639 ± 0.035 
4.591 ±0.059 
4.547 ±0.040 
4.547 ±0.020 
4.740 ± 0.020 
4.555 ± 0.027 
4.567 ±0.021 
4.574 ± 0.026 
4.582 ± 0.033 
4.573 ± 0.028 
4.538 ± 0.045 
4.569 ± 0.042 
4.517 ±0.015 
4.642 ± 0.023 
4.633 ±0.012 
4.610 ± 0.024 
4.629 ± 0.025 
4.721 ±0.047 
4.837 ±0.014 
4.669 ± 0.025 
4.669 ± 0.032 
4.676 ± 0.025 
4.601 ± 0.044 
4.632 ± 0.044 
4.613 ± 0.032 
4.646 ± 0.036 
4.667 ±0.027 



0.006 ± 0.002 
0.009 ± 0.002 
0.008 ± 0.002 
0.007 ± 0.002 
0.009 ± 0.002 
0.007 ± 0.002 
0.007 ± 0.002 
0.008 ± 0.002 

-0.003 ± 0.002 
0.004 ± 0.002 
0.007 ± 0.002 
0.008 ± 0.002 
0.008 ± 0.002 
0.006 ± 0.002 
0.006 ± 0.002 
0.008 ± 0.002 
0.007 ± 0.002 
0.002 ± 0.002 
0.003 ± 0.002 
0.018 ± 0.002 
0.008 ± 0.002 

-0.239 ± 0.002 
0.008 ± 0.002 
0.007 ± 0.002 
0.004 ± 0.002 

-0.164 ±0.002 
0.007 ± 0.002 
0.004 ± 0.002 
0.007 ± 0.002 
0.004 ± 0.002 
0.003 ± 0.002 
0.015 ± 0.002 
0.005 ± 0.002 
0.007 ± 0.002 



Continued on next page. 
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TABLE 1 - Continued 



Morgan et al. 



HJD x 2 /dof Comp. A Comp. B Comp. C Comp. D ref. stars Telescope 



3150 


.512 


1 


.97 


2 


.254 


± 


.008 


2 


087 ± 


.008 


3 


.079 


± 


.012 


4 


.692 ± 0. 


.030 


0. 


.007 ± 


.002 


SMARTS 


3153 


.488 


1 


37 


2 


222 


± 


.008 


2 


110 ± 


.008 


3 


115 


± 


.014 


4 


.708 ± 0. 


046 





.004 ± 


.002 


SMARTS 


3161 


.483 





.83 


2 


.248 


± 


.014 


2 


.089 ± 


.013 


3 


.059 


± 


.022 


4 


.601 ± 0. 


.053 


-0 


.002 ± 


.002 


SMARTS 


3164 


.466 


3 


.56 


2 


.231 


± 


.007 


2 


.170 ± 


.007 


3 


064 


± 


.011 


4 


.632 ± 0. 


.027 





.008 ± 


.002 


SMARTS 


3169 


.555 


2 


90 


2 


256 


± 


.008 


2 


139 ± 


.008 


3 


124 


± 


.013 


4 


.609 ± 0. 


029 





008 ± 


.002 


SMARTS 


3172 


.505 


6 


09 


2 


200 


± 


.007 


2 


196 ± 


.007 


3 


089 


± 


.010 


4 


.602 ± 0. 


.027 


o 


006 ± 


.002 


SMARTS 


3178 


.485 


4, 


37 


2 


.258 


± 


Oil 


2 


171 ± 


.012 


3 


.090 


± 


.018 


4 


.790 ± 0. 


049 





.005 ± 


.002 


SMARTS 


3183 


.473 


1 


.43 


2 


220 


± 


.009 


2 


.234 ± 


.009 


3 


195 


± 


.016 


4 


.557 ± 0. 


038 





006 ± 


.002 


SMARTS 


3185 


.513 


1 


.16 


2 


224 


± 


.009 


2 


.237 ± 


.010 


3 


.198 


± 


.017 


4 


.619 ± 0. 


.043 


0. 


.004 ± 


002 


SMARTS 


3189 


.471 


1 


.36 


2 


.250 


± 


.012 


2 


.251 ± 


.013 


3 


.193 


± 


.022 


4 


.591 ± 0. 


.050 





.002 ± 


.002 


SMARTS 


3192 


.472 


2 


.30 


2 


.275 


± 


.010 


2 


.233 ± 


Oil 


3 


.208 


± 


.017 


4 


.637 ± 0. 


.037 





.007 ± 


.002 


SMARTS 


3201 


474 


1 


.26 


2 


281 


± 


.010 


2 


.261 ± 


Oil 


3 


166 


± 


.017 


4 


.540 ± 0. 


036 





006 ± 


.002 


SMARTS 


3210 


.475 


0. 


.53 


2 


.337 


± 


.032 


2 


.344 ± 


.036 


3 


.279 


± 


.071 


4 


.870 ± 0. 


.231 


— 0. 


.030 ± 


.003 


SMARTS 


3344 


.778 


1 


.47 


2 


.595 


± 


.012 


2 


.805 ± 


.016 


3 


.518 


± 


.026 


4 


.829 ± 0. 


.056 





.002 ± 


.002 


SMARTS 


3348 


.807 


2 


.24 


2 


.614 


± 


.008 


2 


.853 ± 


.010 


3 


.512 


± 


.015 


4 


.808 ± 0. 


033 





006 ± 


.002 


SMARTS 


3352 


.809 


2 


04 


2 


.644 


± 


.010 


2 


.857 ± 


.012 


3 


550 


± 


.018 


4 


.776 ± 0. 


.036 





006 ± 


.002 


SMARTS 


3355 


.807 


2 


59 


2 


.663 


± 


.010 


2 


.903 ± 


.012 


3 


.598 


± 


.018 


4 


.791 ± 0. 


.036 





.006 ± 


.002 


SMARTS 


3357 


777 


3 


.82 


2 


679 


± 


.010 


2 


.933 ± 


.012 


3 


568 


± 


.018 


4 


.769 ± 0. 


035 





.007 ± 


.002 


SMARTS 


3359 


774 


2 


17 


2 


730 


± 


Oil 


2 


900 ± 


.013 


3 


580 


± 


.020 


4 


.750 ± 0. 


.037 





.006 ± 


.002 


SMARTS 


3362 


.784 


1 


.27 


2 


.710 


± 


.010 


2 


.923 ± 


.013 


3. 


.649 


± 


.020 


4 


.789 ± 0. 


.042 


0. 


.003 ± 


.002 


SMARTS 


3368 


.796 


1 


.78 


2 


.748 


± 


.012 


2 


.920 ± 


.015 


3 


.582 


± 


.023 


4 


.756 ± 0. 


.052 


0. 


.004 ± 


.002 


SMARTS 


3373 


.731 


0. 


89 


2 


718 


± 


.014 


2 


881 ± 


.018 


3 


.599 


± 


.030 


4 


.908 ± 0. 


068 





003 ± 


,002 


SMARTS 


3375 


.728 


3 


.54 


2 


.795 


± 


.012 


2 


.829 ± 


.013 


3 


.515 


± 


.020 


4 


.695 ± 0. 


.039 





.005 ± 


.002 


SMARTS 


3380 


.750 


1 


.11 


2 


.710 


± 


.011 


2 


.841 ± 


.013 


3 


.588 


± 


.021 


4 


.753 ± 0. 


.043 


o 


.002 ± 


.002 


SMARTS 


3381 


.762 


1 


68 


2 


721 


± 


.010 


2 


.808 ± 


.011 


3 


567 


± 


.017 


4 


.743 ± 0. 


.033 





005 ± 


.002 


SMARTS 


3383 


.722 


2. 


19 


2 


698 


± 


.009 


2 


834 ± 


.010 


3 


565 


± 


.015 


4 


.763 ± 0. 


032 


0. 


.004 ± 


002 


SMARTS 


3385 


.920 


1 


60 


2 


.741 


± 


.007 


2 


.825 ± 


008 


3 


.633 


± 


.012 


5 


.002 ± 0. 


.030 





.007 ± 


.002 


APO 


3395 


708 


0. 


92 


2 


703 


± 


.014 


2 


842 ± 


.016 


3 


649 


± 


.030 


4 


.828 ± 0. 


.073 


-0. 


.001 ± 


.002 


SMARTS 


3398 


.745 


1 


.20 


2 


.683 


± 


.012 


2 


902 ± 


.014 


3 


580 


± 


.023 


4 


.698 ± 0. 


.051 


0. 


.003 ± 


.002 


SMARTS 


3402 


.736 


1 


.98 


2 


.680 


± 


.009 


2 


.886 ± 


.010 


3 


.636 


± 


.016 


4 


.731 ± 0. 


.035 


0. 


.006 ± 


.002 


SMARTS 


3404 


.833 


2 


.90 


2 


.720 


± 


.010 


2 


.884 ± 


.012 


3 


.575 


± 


.017 


4 


.768 ± 0. 


.035 


0. 


.006 ± 


.002 


SMARTS 


3407 


.773 


1 


46 


2 


733 


± 


.011 


2 


.886 ± 


.013 


3 


664 


± 


.020 


4 


.741 ± 0. 


.035 





.005 ± 


.002 


SMARTS 


3409 


.800 


2 


.36 


2 


.724 


± 


.009 


2 


.935 ± 


.010 


3 


.653 


± 


.015 


4 


.760 ± 0. 


.029 





.005 ± 


.002 


SMARTS 


3411 


.826 


2 


.11 


2 


.741 


± 


.010 


2 


.925 ± 


.012 


3 


.637 


± 


.018 


4 


.724 ± 0. 


.033 


0. 


.006 ± 


.002 


SMARTS 


3414 


.809 


3 


83 


2 


.717 


± 


.009 


2 


967 ± 


.011 


3 


624 


± 


.016 


4 


.786 ± 0. 


033 





006 ± 


.002 


SMARTS 


3420 


.807 


1 


.28 


2 


813 


± 


.012 


2 


.926 ± 


.014 


3 


626 


± 


.021 


4 


.712 ± 0. 


.036 


0. 


.005 ± 


.002 


SMARTS 


3422 


.774 


2 


.40 


2 


.750 


± 


.010 


3 


.021 ± 


.013 


3 


.730 


± 


.020 


4 


.844 ± 0. 


.044 





.003 ± 


.002 


SMARTS 


3430 


.752 


1 


98 


2 


792 


± 


.012 


3 


025 ± 


.016 


3 


631 


± 


.023 


4 


.845 ± 0. 


051 





.004 ± 


.002 


SMARTS 


3432 


.882 


1 


51 


2 


794 


± 


.015 


2 


.975 ± 


.020 


3 


751 


± 


.032 


4 


.823 ± 0. 


.059 


— 0. 


001 ± 


.002 


SMARTS 


3442 


.747 


4, 


.22 


2 


.831 


± 


.010 


3 


.058 ± 


.012 


3 


.696 


± 


.017 


4 


.815 ± 0. 


.034 





.007 ± 


.002 


SMARTS 


3444 


.725 


3 


.22 


2 


.841 


± 


.009 


3 


.066 ± 


.011 


3 


.751 


± 


.016 


4 


.920 ± 0. 


.034 


0. 


.006 ± 


002 


SMARTS 


3446 


675 


2 


62 


2 


868 


± 


.010 


3 


061 ± 


.012 


3 


723 


± 


017 


4 


.876 ± 0. 


.035 


0. 


006 ± 


002 


SMARTS 


3448 


.712 


2 


.75 


2 


.853 


± 


.009 


3 


.087 ± 


.011 


3 


.753 


± 


.016 


4 


.875 ± 0. 


.034 


0. 


.004 ± 


.002 


SMARTS 


3459 


.649 


2 


.73 


2 


.873 


± 


.012 


3 


.063 ± 


.014 


3 


.681 


± 


,021 


5 


.127 ± 0. 


.060 


0. 


.006 ± 


.002 


SMARTS 


3462 


598 


1 


35 


2 


888 


± 


.012 


3 


049 ± 


.015 


3 


724 


± 


.022 


4 


.940 ± 0. 


.043 





005 ± 


.002 


SMARTS 


3462 


.817 


7 


39 


2 


884 


± 


.009 


3 


.003 ± 


.010 


3 


927 


± 


.017 


5 


.166 ± 0. 


.037 


0. 


.002 ± 


.002 


APO 


3464 


.604 


1 


.17 


2 


.921 


± 


.014 


3 


.024 ± 


.016 


3 


.729 


± 


.024 


4 


.955 ± 0. 


.047 


0. 


.005 ± 


.002 


SMARTS 


3469 


.535 


2 


,34 


2 


925 


± 


.010 


3 


083 ± 


.012 


3 


740 


± 


.018 


4 


.975 ± 0. 


.039 





006 ± 


.002 


SMARTS 


3472 


.607 


1 


29 


2 


.970 


± 


.015 


2 


986 ± 


.017 


3 


645 


± 


.024 


4 


.894 ± 0. 


050 


0. 


006 ± 


.002 


SMARTS 


3474 


.744 


5 


,50 


2 


.904 


± 


.012 


3 


.000 ± 


.014 


3 


.714 


± 


.021 


4 


.991 ± 0. 


.047 


0. 


.006 ± 


.002 


SMARTS 


3476 


,683 


2 


,32 


2 


.932 


± 


.011 


2 


.968 ± 


.012 


3 


.699 


± 


.018 


4 


.902 ± 0. 


.038 





.004 ± 


.002 


SMARTS 


3478 


.591 


2 


,55 


2 


879 


± 


.010 


3 


011 ± 


.012 


3 


715 


± 


.018 


4 


.893 ± 0. 


.041 


0. 


.004 ± 


.002 


SMARTS 


3486 


.577 


0. 


,68 


2 


.867 


± 


.022 


3 


.006 ± 


.028 


3 


.680 


± 


.044 


5 


.166 ± 0. 


.121 


-0. 


.008 ± 


.003 


SMARTS 


3491 


.630 


2 


,38 


2 


.865 


± 


.010 


2 


.957 ± 


.012 


3 


.714 


± 


.018 


4 


.927 ± 0. 


.039 


0. 


.004 ± 


.002 


SMARTS 


3501 


.538 


2 


.77 


2 


837 


± 


.012 


2 


.977 ± 


.014 


3 


652 


± 


.020 


5 


.010 ± 0. 


.051 


0. 


.003 ± 


.002 


SMARTS 


3503 


.538 


1 


.20 


2 


.818 


± 


.013 


3 


.014 ± 


.015 


3 


.683 


± 


.023 


4 


.977 ± 0. 


.057 


-0 


.001 ± 


.002 


SMARTS 


3506 


,548 


1 


.81 


2 


.852 


± 


.013 


2 


.983 ± 


.016 


3 


.723 


± 


.026 


5 


.118 ± 0. 


.065 





.005 ± 


.002 


SMARTS 


3508 


521 


0. 


97 


2 


841 


± 


.017 


3 


.001 ± 


.021 


3 


714 


± 


.035 


4 


.922 ± 0. 


080 


-0. 


.001 ± 


.002 


SMARTS 


3511 


,546 


1 


13 


2, 


820 


± 


.013 


3 


.059 ± 


.016 


3 


820 


± 


028 


5 


.075 ± 0. 


070 





001 ± 


.002 


SMARTS 


3512 


,519 


1 


.13 


2 


.824 


± 


.013 


3 


.033 ± 


.015 


3 


.759 


± 


.026 


5 


,165 ± 0. 


.078 





.001 ± 


.002 


SMARTS 


3520 


.601 


1 


18 


2 


930 


± 


.015 


2 


997 ± 


.018 


3 


.725 


± 


.027 


5 


.087 ± 0. 


.061 





.004 ± 


.002 


SMARTS 


3521 


.573 


4. 


54 


2 


861 


± 


.012 


3 


.105 ± 


.014 


3 


.722 


± 


.020 


5 


.046 ± 0. 


.054 


— 0. 


000 ± 


.002 


SMARTS 


3523 


.545 


1 


.70 


2 


.929 


± 


.011 


3 


.036 ± 


.013 


3 


.758 


± 


.020 


5. 


.017 ± 0. 


.043 


0. 


.004 ± 


.002 


SMARTS 


3525 


.568 


2 


.52 


2 


.913 


±0 


.012 


2 


.993 ± 


.014 


3 


.712 


±0 


.021 


4 


.996 ± 0. 


.046 





.006 ± 


.002 


SMARTS 


3528 


.508 


4. 


02 


2 


.893 


±0 


.010 


3 


050 ± 


.011 


3 


752 


±0 


.016 


5 


.079 ± 0. 


042 


0. 


005 ± 


,002 


SMARTS 


3529 


.511 


2. 


.55 


2 


.911 


±0 


.011 


3 


.012 ± 


.012 


3 


.729 


±0 


.018 


5 


.087 ±0. 


.045 


0. 


.005 ± 


.002 


SMARTS 


3542 


.517 


1 


.00 


2 


.855 


±0 


.024 


2 


.837 ± 


.026 


3 


.653 


±0 


.047 


5 


.104 ±0. 


.137 


-0. 


.012 ± 


,003 


SMARTS 


3550 


.558 


0. 


60 


2 


810 


±0 


.027 


2 


672 ± 


.026 


3 


.575 


±0 


.048 


5 


.202 ± 0. 


155 


-0. 


019 ± 


,003 


SMARTS 


3569 


.458 


1 


25 


2 


599 


±0 


.013 


2 


602 ± 


.014 


3 


.468 


±0 


.025 


5 


.268 ± 0. 


.107 


0. 


001 ± 


,002 


SMARTS 


3570 


.493 


1 


.07 


2 


.623 


±0 


.022 


2 


.531 ± 


.022 


3 


.464 


±0 


.042 


5 


.361 ±0. 


.171 


-0. 


.012 ± 


.003 


SMARTS 


3578 


.475 


1 


58 


2 


.590 


±0 


Oil 


2 


521 ± 


.012 


3 


.431 


±0 


.019 


5 


.040 ± 0. 


055 


0. 


004 ± 


.002 


SMARTS 


3580 


.478 


1 


.28 


2 


.572 


±0 


.019 


2 


527 ± 


.020 


3 


395 


±0 


,036 


5 


.854 ±0. 


238 


-0. 


013 ± 


.003 


SMARTS 


3582 


.478 


1 


.87 


2 


.607 


±0 


.026 


2 


.458 ± 


.025 


3 


.352 


±0 


.040 


5 


.370 ± 0. 


.140 


-0. 


.004 ± 


.002 


SMARTS 


3583 


.461 


1 


.12 


2 


.558 


±0 


.024 


2. 


.500 ± 


.026 


3. 


.321 


±0 


.044 


5 


.549 ± 0. 


.220 


-0. 


.014 ±0 


.003 


SMARTS 



Note. — HJD is the Heliocentric JD minus 2450000. The A-D columns give the image magnitudes relative to a local 
reference star. The ref. stars column gives the average magnitude differences relative to the campaign mean. 
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TABLE 2 

HST ASTROMETRY AND PHOTOMETRY OF RX J1131 — 1231 



Comp. 


ARA 


ADec 


V 




H 


A 


= 


= 


18.20 + 0.14 


17.77 ±0.10 


15.97 ±0.02 


B 


+0.031 ± 0.003 


+1.187 ±0.003 


18.25 ±0.03 


18.11 ±0.11 


16.22 ±0.04 


C 


-0.588 ±0.001 


-1.120 ±0.000 


18.51 ±0.05 


18.19 ±0.02 


17.02 ±0.02 


D 


-3.105 ±0.003 


+0.879 ± 0.002 


20.02 ±0.12 


19.81 ±0.03 


18.76 ±0.03 


E 


-1.933 ±0.003 


+1.139 + 0.004 


24.47 ±0.22 


22.73 ±0.02 


20.88 ±0.08 


G d 


-2.032 ±0.002 


+0.586 ± 0.001 


19.63 ±0.06 


17.88 ±0.03 


16.13 ±0.24 


Go 


= -2.032 


ee +0.586 


21.72 ±0.53 


19.39 ±0.35 


17.29 ±0.14 



Note. — Relative astrometry and absolute photometry of RX J1131— 1231 com- 
ponents from HST/ACS (V- and /-band) and ffST/NICMOS (ff-band) observations. 
Relative positions arc from the H-band data. Components and G e denote the dc 
Vaucoulcurs and Exponential profiles for the galaxy model. 



TABLE 3 
Model Results for RX J1131-1231 



Model 


N do f 


x 2 


Xpos 






b 


G x 


Gy 


e 


9e 


7 


6-y 


s 


TA-B 


TA-C 


TA-D 










(") 


(") 


(") 










(") 


(days) 


(days) 


(days) 


SISx 


6 


538 


475 




63.4 


1.86 


EE2.032 


=0.586 






0.155 


-73.6 




0.86 


0.99 


-119 


SIEx 


4 


138 


77.2 




60.8 


1.82 


=.2.032 


=0.586 


0.182 


-57.0 


0.112 


-84.9 




0.98 


1.25 


-117 


SIEx+ 


4 


74.1 


2.9 


10.1 


61.0 


1.83 


2.036 


0.571 


0.162 


-59.1 


0.113 


-82.6 




0.98 


1.23 


-117 


SIEx+ 


2 


23.8 


0.1 


0.2 


23.5 


1.55 


2.033 


0.585 


0.132 


-62.6 


0.035 


81.7 




5.90 


7.88 


-80 


a = 1 












0.32 


-0.049 


-0.035 


= 








0.16 








SIEx+ 





4.7 


0.2 


0.2 


4.3 


1.66 


2.032 


0.586 


0.165 


-63.4 


0.052 


87.4 




12.00 


13.78 


-85 


a = 0.72 












0.27 


-0.056 


-0.041 


= 








= 0.2 









Note. — Best-fit model parameters for five models described in the text. For the a — 1 and a. — 0.72 perturber models, the first line gives the 
model parameters for the primary SIEx+ halo and the second line gives the model parameters for the perturber near image A. All relative positions 
are measured with respect to image A. 



